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This paper is concerned with phase contrast in electron images with 
emphasis on a periodic scattering object. The Kirchoff diffraction or imaging 
integral over the back focal plane is formulated in terms of the amplitude and 
phases of the scattered wave, using coordinates adapted to numerical methods. 
The integral was programmed for evaluation on the IBM 7090 and the 
image plane amplitude displayed on a microfilm plotter. The effects of 
spherical aberration, aperture size, defocus and thermal motion of the scatter- 
ing atoms on image contrast of atom positions for chains of nickel and of 
gold atoms were investigated for separations of 2 A and 8 A. 

In the range of atomic separations, spherical aberration is the most de- 
structive single factor in the loss of phase contrast. It appears that an ob- 
jective lens with a spherical aberration coefficient less than 0.2 mm will be 
necessary if phase contrast images of atom locations are to be attained. In 
addition, a practical quarter-wave phase plate is essential for the objective 
lens system and will be much more effective than defocus contrast. Even so 
the contrast is marginal for photographic recording except for the case of a 
thin perfect crystal. Amplitude eotitrast is very small for atom positions and 
should be minimized by (he use of a large objective aperture. 

Phase contrast shoidd improve with increasing accelerating potentials due 
to reduced inelastic cross sections compensating the loss in elastic contrast 
near 2 X W 5 volts. 
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I. INTRODUCTION 

The resolving power and general quality of electron microscope images 
are determined basically by image plane contrast. There are two distinct 
contrast mechanisms: 

(1 ) amplitude contrast — produced by removal from the image plane 
of electrons scattered outside the objective half -angle, /3 bj 

(2) phase contrast — produced by suitable recombination at the 
image plane of waves scattered within the objective aperture. 

The first mechanism, amplitude contrast, is the one commonly operat- 
ing in images of objects which are greater than about 10 A in size. The 
contrast between two image points for this type of object is either mass 
thickness or diffraction and is given approximately by 

AI/I = A(Qt) = QAt or tAQ. (la) 

for amorphous materials. (See Ref. 1, Ch. IX, for diffraction contrast.) 
Here Q is the cross section for scattering outside the objective aperture 
and t is the object thickness. In turn, the cross section for scattering out- 
side the aperture is 

Q = AW™ p (lb) 

with No = 6.02 X 10 23 being Avogadro's number. Here, A is the atomic 
weight and p the density. The cross section per atom for scattering out- 
side the aperture is o- atom and consists of an elastic and an inelastic part 

= O-el + 0"inol • (1C) 



Catc 



The cross section (r atom is the fraction of incident electrons scattered 
beyond the objective half-angle, /3 bj , and is the integral of the differen- 
tial cross section, Z>(/3), over the scattering angle /3 

<r„tom = 2t [ Z>(/3) sin <#. (Id) 

The differential cross section, Z>(/3), determines the intensity distri- 
bution at the back focal plane of the objective lens. For elastic scattering 
by an isolated atom D(/3) is simply the square of the atomic scattering 
amplitude or 

D(0) = \f(s)\ 2 (2a) 

where /(s) is the electron scattering amplitude per atom as a function of 
the scattering parameter 

s=4t! M0/2) 
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with the scattering angle and X the incident electron wavelength. For 
the small scattering angles encountered with fast electrons, sin (0/2) « 
0/2 and 

.«4,|. (20 

Tables of scattering amplitudes [Ref. 2, Tables 3.3.3 A(l) and A (2)] 
list values of f(s) in A or cm for chosen values of the parameter 0/2X. 

The differential cross section for a solid object is determined by inter- 
ference among the waves scattered by the individual atoms composing 
the object. This, in turn, depends upon the spatial arrangement of the 
atoms. The differential cross section for a liquid or glass is diffuse, with 
broad maxima corresponding to those interatomic distances occurring 
most frequently. On the other hand, that for a single crystal is a discrete 
set of spots. 

The portion of the scattered amplitude distribution falling outside 
the objective aperture determines the amplitude contrast of an image 
point as sketched in Fig. 1. That part of the distribution falling within 
the aperture is available for phase contrast at the image plane. The 
extent to which the phase information can be used in the formation of 
an image is presently limited by spherical aberration which essentially 
"scrambles" the phases to a degree increasing as the fourth power of the 
scattering angle. Phase information concerning the distance between two 
object points is thus increasingly garbled as the distance between the 
points decreases. This follows since the scattering angle for a distance 
a in the object is 

« X/o. (2d) 

Only values of a such that < o i,j have any possibility of contribut- 
ing to phase contrast. 

Although numerical calculation of phase contrast is the purpose of 
this paper, it is instructive to consider the limitations on amplitude 
contrast in the region of interatomic separations before discarding it. 

II. AMPLITUDE CONTRAST 

The total amplitude contrast for a single atom in the object plane of 
a perfect lens would be 

— ~ - (r "' _ (rinp ' (3 a ) 

Zo ~ tt(RJ) t(RJ) 
where R v \ is the scattering radius for clastic collisions and fi\. x that for 
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Fig. 1 — Schematic diagram showing amplitude scattered by single atom in 
object plane. The diffracted amplitude at the back focal plane produces image 
plane contrast depending upon the size of the aperture. The portion (shaded) 
falling outside the aperture results in deficiency amplitude contrast. The portion 
within the aperture can produce phase contrast. 
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inelastic. If the resolving power of an imperfect objective lens is 8, then 
R e i £» S in (3a). Both terms of (3a) rapidly diminish with increasing 
incident electron velocity and fixed aperture. The negative sign denotes 
deficiency contrast for an image point. 

Since the scattering contracts into the forward direction as the 
electron velocity increases, it is advisable to consider the scattering 
amplitude at constant values of /3 hj/2\ (the scattering parameter). 
This amounts to shrinking the objective aperture as the wavelength de- 
creases so that the Airy diffraction disc from the aperture maintains a 
constant radius. This is accomplished by changing the variable of inte- 
gration in (Id) from /3 to the scattering parameter s of (2b). The elastic 
cross section for scattering outside the objective aperture (per atom) 
becomes 

,2/, 2i-l fir/X 

a Bl = X(1 7 T) / \f(s)\ 2 sds (3b) 

with 7 = v/c, the ratio of electron velocity to that of light. Here / (s) 
is the scattering amplitude based on the electron rest mass. It is noted 
that the velocity dependence of (3b) results in a rapidly decreasing value 
of o-ei as the accelerating voltage increases, with consequent loss of elastic 
amplitude contrast in (3a). The elastic cross sections as a function of 
atomic number have recently been calculated by Burge and Smith 
using numerical methods with (3b ) . 

The second term of (3a) depends upon the inelastic cross section 
(Ti ne i and the inelastic excitation radius, R PX . Williams 4 arrived at an 
excitation distance R ex for an atom of average excitation energy (AE) 
having a velocity dependence 

0.2vh , 2,-\ , ,s 

R ~~W) ( y) ■ () 

For carbon, R cx is about 40 A at V„ = 10 5 volts taking (AE) tt 44 volts, 
so that the location of the energy loss event is very diffuse. For heavier 
atoms, .Rex is generally less, being about 4 A for Al at 10 5 volts with 
(AE) ^ 300 volts. The inelastically scattered electrons thus do not carry 
localized information on atom locations (unless they are subsequently 
diffracted). The relatively large values of R ex make the amplitude con- 
trast due to inelastic scattering trivial except for relatively large objects 
(50-100 A). 

The amplitude contrast resulting from elastic scattering is so small as 
to offer little hope of imaging single atoms. For example, the elastic cross 
section for a nickel atom is about 7 X 10 18 cm 2 at V a = 5 X 10 4 volts 
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for an objective aperture /3 bj «5X 10~ 3 rad. The contrast is only about 
0.8 per cent for 8 tt 2 A, or so low as to be of little importance. The 
chief hope for imaging atoms with the electron microscope lies in phase 
contrast, wherein intensities are very sensitive to the relative phases of 
the scattered waves reconstituted at the image plane. 

III. PHASE CONTRAST 

The dominant contrast mechanism for scattering objects exhibiting a 
periodic structure or with a detail size approaching atomic dimensions 
is that of phase contrast. This requires that the image plane amplitude 
be evaluated by superposition of all the waves which are scattered by 
the object and enter the objective aperture of the lens. In its most general 
form this superposition is the imaging integral. 5 For a perfect, aberration- 
free objective lens the imaging integral is the Fourier transform of the 
diffracted amplitude distribution at the back focal plane, which is itself 
the Fourier transform of the distribution of scattering potential in the 
object. The imaging integral is thus a magnified representation of the 
Fourier projection in the object plane of the electrostatic potential distri- 
bution 6 in the object itself. In the case of a real objective lens, the 
fidelity of the image plane amplitude distribution is determined by the 
lens aberrations and by the size of the objective aperture. The objective 
aperture limits the distance information available to the image plane, 
removing the higher "spatial frequencies" in the diffraction pattern at 
the back focal plane. This is in addition to its own diffraction pattern. 

A description of the intensity distribution in the image plane involves 
a consideration of three planes in the objective lens. These are: (1) the 
object plane with rectangular coordinates, the x°-y° plane, (2) the back 
focal plane or £'-77' plane containing the objective aperture and (3) the 
image or x'-y' plane. The three planes are normal to the optic axis z 
as shown in Fig. 2. The object and image distances measured from the 
lens plane are L and L, respectively. The distance from the lens plane 
to the back focal plane is the focal length of the lens /. The magnifi- 
cation of the image is M = Li/L , and / tt Lo when L, ^> L . 

Assume a plane parallel electron beam of wave vector K along the 
optic axis incident on the object plane from above. Let its amplitude be 
unity. Electron waves scattered by the object in the object plane pro- 
duce a diffraction or scattering pattern at the back focal plane. This 
pattern is nearly a plane section through the reciprocal lattice of the 
object with a scale factor L X, where A is the electron wavelength and 
| K | = 2rr/\. The diffracted amplitude distribution ^ in the back focal 
plane is in turn a source for Huygens wavelets which propagate to the 
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Fig. 2 — The three planes — object, back focal and image — that are of con- 
cern in contrast calculations. The scale in the back focal plane is for a periodic 
scattering object with reciprocal lattice vectors g\ and g» . 

image plane and recombine according to their phases and so produce the 
image amplitude distribution. The entire process can be formally de- 
scribed in terms of Fourier transforms, as already mentioned. However, 
when lens aberrations are present the imaging integral describing the 
image plane amplitude distribution for the imperfect image is no longer 
identical with the object plane Fourier projection of potential in the 
object. It is usually necessary to employ numerical methods in evaluating 



214 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1965 

the image plane amplitude, since the integration cannot in general be 
carried out in analytic form. The integral to be considered 7 is 

*(xW) = A^ ( f *<£*> ex P MM)] 

•exp [-iKiit'x' -t.vV)/^] di dr,'. 

x($', v) is the phase of the amplitude at ((•', ?/) in the back focal plane. 
Since the primary concern is with crystalline or periodic scattering 
objects, the coordinate system chosen in the back focal plane is adapted 
both to this situation and to a formulation lending itself to numerical 
evaluation. Distances R in the back focal plane are related to distances 
d in the object by the relation | R | = L X/d. L \ is the camera constant 
or scale factor. A single-crystal object produces an array of diffraction 
spots in the back focal plane located at R = L \gi + L \g2 where gi 
and g 2 are the operating reciprocal vectors in the object plane. The 
generating vectors in the object are ai ,a 2 ,a 3 and the generating vectors 
in reciprocal space are bi ,b 2 ,b 3 , subject to the condition 

(a,- by) = 8ij 

= if %9*j (6a) 

= 1 if i = } 

so that | bi | = | a x | — etc. and 

g = M>i + fcb 2 + Zb 3 (6b) 

where h,k,l are the Miller indices denoting the reciprocal lattice vector g. 
Let £ and 77 be numbers such that an arbitrary point in the back focal 
plane is the terminus of the vector R 

R($', 7/) - £LoA gl + nLoAfc . (6c) 

If £ and 77 are integers, R locates a diffraction spot at a distance 

|fi|- [(&#*? + (r/LoA^) 2 ]* (6d) 

since the £ and 77 axes are rectangular. The scattering angle fi for electrons 

to the point R is, for small angles (sin /3 « /3) 

- 1*1 - xK&i) 1 + (n*)V (6e) 

and the scattering parameter s is 

s . *I S in j? « Mi&x)* + (W)*?- (Of) 

A ^ 
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The amplitude at a point (a;*, y l ) in the image plane is the diffraction 
integral 8 over the aperture 

/•+max 

L L,X 

•exp [ — 2iri(£c o 0i + WWl d(Lo\£gi)d(L \r}g2) 

LoL.X ■'-u.x-'-^., (7 a ) 

•exp [— 2*t(£r°0i + flfa) d£ di?. 

In order to evaluate (7a) it is necessary to know the amplitude distribu- 
tion iH£, v) and the phase x(£, v) in the back focal plane. 

If the diffraction pattern has a center of symmetry at the origin of 
the back focal plane, the periodic object likewise has a center of sym- 
metry at the origin of the object plane. In this case the imaging integral 
(7a) can be simplified for the purposes of calculation to 

*(*V) = ^ /*"(*'"*(**)[«» x (fc») + I'dnxtt*)] 

M J o J o 
X cos 2-tr^xgi cos 2mry°gd£ d-q. (7b) 

_ 2X0102 ^i 
M 

where »S" denotes the phase integral in (7b). 

It is now necessary to choose expressions for the amplitude ^(£, rj) at 
the back focal plane occurring in (7b). This amplitude depends upon 
the detailed spatial arrangement of atoms in the object plane and upon 
the atomic scattering amplitude f(s). There are two approaches to 
^(£, v) — the kinematic and the dynamical theories of electron dif- 
fraction. 9 

The kinematic approach is much the simpler of the two and will 
illustrate the pertinent features in phase contrast. Since the main 
interest here is with scattering objects which are basically periodic, the 
kinematic amplitudes at the back focal plane can be immediately written 
down. Let the object be a crystal sheet with an atom at the origin of 
coordinates in the object plane. If the lateral extent of the crystal is wi 
atoms along x° and n 2 atoms along if and if the extent in the z direction 
is r? 3 atoms, the kinematic amplitude is the familiar expression 10 

./„ v ( , 2x-MV \ sin tttii^ sin im-m fa \ 

MS,n) = (1 - 7 ) T(s)n 3 —. : (8) 

T vs " " sin irf sin ttti 



sin 7r£ sin 7tt? 
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where Wi , n 2 , and n 3 are odd integers. If even integers are employed, 
phase factors must be included in (8). The discussion is therefore con- 
fined to odd values. The atomic scattering amplitude / (s) is based on 
the rest mass of the electron and (1 — y 2 )~ is the relativistic correction 
(with 7 the ratio of electron velocity to that of light). It is noted that 



K 



and 



7Tttl£ = !~jr-! 7ll — 

^ 01 



\K\ A 

irntq = I-J5-! ri2 — . 



2 <7 2 
The phase x in the back focal plane is made up of several terms 11 

The phase change upon diffraction by the object is 7r/2, while the second 
term is the phase due to spherical aberration. The third term is the phase 
introduced by defocusing the lens an amount A/. There is a fourth term 
which is the phase associated with the electron scattering amplitude /(s) 
for an atom, but it is small except for the heaviest atoms and at large 
scattering angles and will be neglected here. 

It should be mentioned that the spherical aberration phase term in 
(9a) is four times larger than that commonly used in discussions 1213 of 
the effect of spherical aberration. The term in (9a) applies to a single 
ray at angle /3 with a corresponding circle of aberration of radius C a /3 . 
The unweighted average radius for a bundle of rays filling the range 
< < /W is 



3 
max 



<a/). p = S±.f aix f?dfi-\c# 

/3 m ax J ° 4 

with a resultant average phase 

(x). p = y^T^- (9b) 

The use of (9b) in (9a) in evaluating the imaging integral is not correct, 
since it suppresses the destructive effect of spherical aberration. 

The total amplitude at an image point is the sum of the scattered 
amplitude (7b) and the unscattered axial wave. If the incident amplitude 
is unity, then in the kinematic approximation of weak scattering the 
unscattered wave leaving the object plane is very nearly of unit ampli- 
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tude. Its amplitude at the origin of the back focal plane is large, but at 
the image plane it is (M)~ l . Adding the unscattered amplitude to (7b) 
yields the total image point amplitude 



*(*', »')«.! = li [l + 2X^2(1 - Y 2 r*£']. 



(10a; 



The mass correction factor (1 — 7 ) ~' has been removed from S l so that 
the velocity dependence of the amplitude can be easily seen. The expres- 
sion (10a) does not apply when dynamical conditions 14 are realized in 
the object, since then the axial wave may be weaker than the scattered 
amplitude and the approximation is no longer valid. 

The image point intensity obtained by multiplying (10a) by its 
complex conjugate is 



VixWfl 



M 2 



1 +4X0102(1 - 7 2 ) te'reall 



iob: 



since the term in | S* \ 2 will generally be trivial compared to the cross 
product term until the mass correction term becomes large at very high 
accelerating voltages (£^10 6 volts). The imaginary part of the scattered 
amplitude (7b) is here neglected. 

The kinematic contrast G between two image points is defined to be 
the intensity difference between the respective points divided by the 
background intensity. Using (10b) the result is 

G&4tafe(l -tVaS* (10c) 

with all lengths in angstroms. A.S" is the amplitude differential between 
the image points in question obtained by numerical evaluation of 

Mm.1 rwx sin rii7r£ sin mm 

& = n 3 / / / (s) — : 

Jo J o sin tt| sin m (lOd) 

X cos x (£,*?) cos 2ir£cV cos 2niy O 0« d£ drj 

at each point. 

For the purpose of phase contrast calculations it is sufficient to con- 
sider a single chain of atoms of spacing a lying along the x axis. The 
number of atoms in the chain is n\ , with the middle atom on the optic 
axis. Setting n 2 = n 3 = 1 in (lOd) gives the integral to be evaluated at 
a point on the x' axis as 

Jo J o sin 7r£ \a / 
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having maxima at the image points corresponding to atoms or where 
x° = 0, ±o, ±2a, • • • etc. The amplitude maxima for atomic positions 
approach 

r«-„ p« / 0( s ) ^1^1 cos x(fo)d* dr, -> rFf(s) cos x(f,*)d* drj 

Jo J o sin x? ■'o •'o 

as the size of the objective aperture, £ mtt x , imax , increases. The value of 
the maxima is seen to be sensitive to the phase x(£, v)- For a perfect 
lens, the scattered amplitude in the Gaussian image plane vanishes, 
since x ~ t/2, so there will be only background intensity. On the other 
hand, the introduction of a quarter-wave phase plate 15 into the scattered 
beams can make x = or ir and the amplitude will be a maximum. The 
sign of the maxima will be that of cos x, so that atom positions can be 
either bright or dark against the background, depending upon the phase. 

IV. NUMERICAL RESULTS 

4.1 Spherical Aberration at the Gaussian Image Plane 

Numerical evaluation of the phase integral S % in (7b) requires that 
the phase (9a) be expressed in terms of the dimensionless coordinates 
(£, 7] ) using (6e ) to give 

xfc, V) = (t/2) - 27rC 6 \ 3 [(& 1 ) 2 + (m?f + ttA/XK&O 2 + (Mt)*\. 

For the single-atom chain of spacing a this reduces to 

xft, v) ~ (t/2) - *C 6 '[? + v? + rAf(f + t, 2 ) (11a) 

with C 6 ' = 2C X 3 /« 4 and A/' = A/A/a 2 being dimensionless aberration 
and defocus parameters. It is noted that %C d ' is the number of wave- 
lengths of spherical aberration at the first diffraction maximum (£ = 1, 
■q = 0) while |A/' is the number of wavelengths of defocus. 
In the Gaussian image plane A/ = and (11a) becomes 

X&v) = (x/2) - irC 6 '[l? + v 2 ] 2 . (Hb) 

The real phase integral <S* along the chain, y° = 0, is now 

S i = f "■« /•'" ax ^^ cog x( ^ ) cog 2v $ x «/ a dk dn (12) 

J •'0 

and the absolute amplitude is 2\a~ 2 S\ If the phase (lib) is used in 
(12) it is evident that cos x = sin wC 6 '[£ + rj 2 ] 2 . The trivial contrast 
observed is produced by spherical aberration. Fidelity contrast at the 
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Gaussian image plane requires the use of a A/4 wave plate in the back 
focal plane. If the plate advances only the diffraction spectra but not the 
unscattered or zero-order beam, the phase term in (12) is 

cos x = — cosirGoUf + V J 

and the image will display atom positions dark relative to the back- 
ground. 

Since analytical expressions for the atom scattering factors are not 
available, the amplitude in the back focal plane for this purpose is ob- 
tained numerically from tabular values as described in the Appendix. 
A discussion of the choice of sampling intervals along the £, tj and x /a 
axes is also found in the Appendix and some of the artifacts that may 
occur are pointed out. 

If Cq = and n x = 1 in (12), it will be noted that S l is the profile 
of the Airy disk 16 for a single atom in the object plane due to diffraction 
by the objective aperture. Fig. 9(b) of the Appendix shows a computer 
plot (IBM-7090 microfilm plotter) of (12) for this case of a single 
nickel atom. The absolute maximum amplitude is 7.7 A angstrom, or a 
contrast of about 60 per cent relative to background for 100-kv elec- 
trons. The value a = 2 A was employed in (10c), since the scaled /°(s) 
curve was at this value of a. The contrast is independent of the number 
of atoms, ni , in the chain but is proportional to n A , the number of atoms 
along the optic axis if chains are stacked one directly above the other. 

Since the electron microscope objective aperture imposes a finite 
upper bound on the imaging integral the result is equivalent to termi- 
nating or truncating a Fourier series to a finite number of terms. As a 
consequence, the imaging integral does not converge to the value ex- 
pected for an infinite upper bound but oscillates about this value. The 
Gibbs phenomenon 17 in the particular case of the objective aperture 
manifests itself as the experimentally observed diffraction pattern of the 
aperture. The same phenomenon gives rise to artifacts between atom 
positions in the image when a finite number of diffraction spectra are 
admitted by the aperture. 1 

The devastating effect of spherical aberration upon phase contrast in 
the image of an atom chain is illustrated in the series of amplitude pro- 
files of Fig. 3. The numerical values are for nickel with a spacing of 2A. 
Profiles were obtained for chain lengths of 5 and 17 atoms which are 
within the transverse coherence 19 length of double condenser illumina- 
tion in present microscopes. Inelastic scattering and thermal motion are 
neglected. The former is reduced by increasing the accelerating potential 
and the latter by reducing the temperature. The profiles of Fig. 3 speak 
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Fig. 3 — Series of numerical image amplitude profiles for a 5-atom chain of 
nickel atoms with spacing a = 2 A. The effect of increasing spherical aberration 
parameter Cq' (£ Cq' is the number of wavelengths of aberration at the first diffrac- 
tion maximum) is illustrated in the series (a)-(e): (a) Cq = 0, (b) 0.2, (c) 0.3, 
(d) 0.5, and (e) 1.0. A quarter- wave phase plate is assumed present in the back 
focal plane. The contrast of about 60 per cent in (a) falls rapidly to about 9 
per cent in (b) for 100-kv electrons with/ = 3 mm. 
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for themselves, showing a loss in atom position amplitude from 15.5 A 
with C 6 = to 2 A with the introduction of 0.1 wavelength of spherical 
aberration at the first diffraction maximum where C 6 = 0.2. 

If the diffraction spectra in the back focal plane were points, as for an 
infinite, perfect crystal, the phase contrast would be periodic with C 6 ' 
and would have maxima where (Az°) 8p h = Cg/3 3 is an integral number 
of half-lattice spacings. For the finite chains of 5 and 17 atoms, the 
damped periodic nature of the contrast is illustrated in Figs. 4 (a) and (b) 
compared to that for a single atom. The aperture has been decreased 
from 3.5 diffraction maxima in Fig. 4(a) to 1.5 in Fig. 4(b) with a 
sizable reduction in contrast. The reduction in aperture size further 
limits the amount of spherical aberration but this is more than compen- 
sated for by the loss of distance information in the diffraction pattern 
and the increased diffraction by the aperture. It must be concluded that 
if useful phase contrast of atom positions is to be obtained spherical 
aberration must be minimized. Other contrast enhancing devices cannot 
overcome this defect of disturbed phase information. 

4.2 Contrast by Defocus with Spherical Aberration 

A case approximating the present state of the electron microscope ob- 
jective lens is that for no phase plate with C 6 ~ 3 mm and contrast en- 
hancement by defocus. This example is a chain of 5 gold atoms with a 
spacing of 8 A and was chosen to approximate the situation of gold- 
stained sites on a DNA molecule. 20 No allowance is made for a substrate 
or inelastic scattering. 

The phase (11a) for this particular case retains tr/2 (with no phase 
plate) as well as the defocus term. The scaled aberration parameter is 
now C 6 ' = 0.74 or 0.37 wavelengths at the first diffraction maximum 
using 100-kv electrons. A series of amplitude profiles was computed on 
the IBM 7090 for a range of values ^ A/ ^4. From these results the 
relative amplitude of atom positions A*S" was plotted against Af as 
shown in Fig. 5. The oscillation of contrast with changing defocus is 
typical with near zero contrast at exact focus. Maximum contrast is ob- 
tained by weakening the objective lens (Af positive) as is well known. 
Even so the contrast does not rise above about 6 per cent, which is sub- 
marginal for seeing the gold atom positions in the image. 

If the spherical aberration coefficient is reduced to 2 mm or C 6 « 0.5 
with 0.25 wavelengths at the first diffraction maxima, the maximum de- 
focus contrast rises to around 9 per cent for a gold atom with 100-kv 
electrons. The neglect of substrate, thermal motion and inelastic scatter- 
ing again renders the visibility of single gold atoms marginal at best. 
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Fig. 4 — Plot of AS' vs Co' = 2C 6 \ 3 /a* for nickel atom chains of 5 and 17 atoms. 
The curves approximate damped, periodic functions. The kinematic contrast is 
G = 4.4 AS' per cent for / = 3 mm using 100-kv electrons. The upper bound 
is 3.5 diffraction maxima in (a) and 1.5 in (b). 
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Fig. 5 — Defocus phase contrast for a chain of 5 gold atoms with 8 A spacing 
and a spherical aberration parameter ('§' = 0.74 or 0.37 wavelengths at the first 
diffraction maximum (Ca = 3 mm). The defocus is {-Af wavelengths or A/ = 1730 
A/' angstrom. The contrast is G « 0.28 AS' per cent at 100 kv. 

However, a clump of three gold atoms at the staining sites might be 
visible in the image, since this would multiply the contrast by a factor of 
roughly two. There might be also some benefit from amplitude contrast 
providing the phase and amplitude contrast are the same sign and that 
the image points for the two coincide. A small shift in image points by 
unsymmetric phase would separate the phase and amplitude image 
points, resulting in confusion. 

If the diffraction spectra at the back focal plane are discrete points, 
the condition for maximum phase contrast by defocus is that the phase 
(11a) be rnr where n is an integer. The optimum defocus parameter A/ 
is then 

A/' = (n - *) + C 6 ' 

a relation useful in estimating the amount of defocus for best contrast. 
Because defocus can optimize only one object plane spacing at a time, 
a phase plate is much to be preferred wherein all spacings are maximized 
in the same image plane. 



4.3 Effect of Thermal Motion 

The scattering amplitude ^(£, -q) in (12) assumes that the atoms in 
the object plane are stationary. This is not true, since they possess ther- 
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inal motion which is temperature-dependent and zero-point motion at 
the absolute zero of temperature. Detailed analysis of thermal vibration 
amplitudes in a solid is a complex problem which need not be discussed 
here. To a first approximation the effect of an isotropic thermal motion 
is described by the Debye- Waller factor e~ M to produce an effective atom 
scattering amplitude 

(1 - tVAsK" 

This dependence has been recently experimentally verified by Horstmann 
and Meyer. 22 The uncertainty lies in the evaluation of M. For the iso- 
tropic averaged vibration case, M is given by 

^ (« 2 > sin 2 0/2 « 2tW>[(^i) S + (mh)*\ 

for small angles. For the simple one-dimensional grating (u a ) is the 
mean square atom displacement along the chain. The Debye factor 
is then 

M«2x 2 ^- > (^ 2 + r 7 2 ) (13) 

a 2 

and the effect of thermal vibrations on contrast can be approximated by 
introducing e~ M into the phase integral (12). The thermal diffuse am- 
plitude is neglected, as is the inelastic scattering, so that again the com- 
puted contrast will be higher than could actually be expected. The effect 
on contrast of thermal motion alone arising through diminution of the 
diffracted amplitudes is illustrated using (12) and setting C Q = A/ = 
and inserting a quarter- wave plate in the back focal plane. Amplitude 
profiles for a range of values of the relative mean square displacement 
were computed. The results are summarized in Fig. 6, showing A*S* as a 
function of (uf)ya. The contrast for stationary atoms is 60 per cent, as 
before, but falls to 30 for a relative mean square thermal displacement of 
0.1. If spherical aberration were introduced, the contrast of atom posi- 
tions would rapidly fall below that necessary for visibility in the usual 
microscope viewing system. 

v. DISCUSSION 

The foregoing numerical results serve to point up the rather severe 
requirements for an objective lens system capable of yielding phase con- 
trast images of atoms. The highly destructive influence of spherical 
aberration is a major hurdle that must be reduced to a minimum. If an 
objective lens with C « 1 mm can be realized, the best contrast obtain- 
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Fig. (i — Effect of thermal motion on image profile maximum amplitude; 
(h,, 2 ) * is the root mean square amplitude along the chain. The values are for 
nickel atoms with a = 2 A and (',-/ = A/' = 0. 



able with an ideal phase plate will be around 10 per cent for an individual 
nickel atom in a chain. The situation is much more favorable for a thin 
crystal of the order of 50 A thick, wherein the diffraction spectra arc 
greatly reinforced. For the single layer of atoms it may be necessary to 
turn to image intensifiers in combination with communication tech- 
niques of extracting a useful signal from background noise. 

Practical considerations demand ultra-high vacuum to eliminate con- 
tamination by electron bombardment. A cryogenic stage will be useful 
to reduce loss of contrast due thermal motion. The problem of back- 
ground noise in the image from a substrate can be circumvented either 
by a self-supporting specimen over a small hole or by using thin single- 
crystal substrates. The noise level from carbon substrates is intolerably 
high for this purpose. 

The rapid decline of inelastic scattering cross sections with increasing 
accelerating voltage should offset the reduced elastic contrast predicted 
by (10c), suggesting that potentials in the range 150-200 kv should be 
appropriate. This is just opposite to the use of lower electron velocities 
for amplitude contrast. 

It appears necessary to develop a practical quarter-wave contrast plate 
for the back focal plane. The thickness t of a region of potential V re- 
quired to introduce a phase advance of x/2 is" 



ttTV 



" ± 

= - or t = 



2V 



where V a is the accelerating voltage. A material film of inner potential 
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V is a possibility, but the elastic and inelastic scattering by the film 
itself must be reckoned with. 

The authors would like to acknowledge the invaluable aid of Miss 
Barbara Dale, who programmed the imaging integral and carried out 
the computational procedures and checks. 

APPENDIX 

The appropriate numerical values of the amplitude in the back focal 
plane are obtained from tabular values of the atom scattering factors. 
The scale factor for/°(s) in the back focal plane requires that 

/ t% , , 0n sin <8/2 
V £ 2 + v 2 = 2a — i — • 

For a nickel atom chain with a = 2A , the first diffraction maximum is 
at (sin/3/2)/X = 0.25 and Vt 2 + V 2 = !> for whic h /" = 3.22 A^corre 
sponding to a scattering angle /3 = X/a V£" + v 2 = 1-8 X 10 rad. 
The empirical curve /"(Vii 2 + V 2 ) ls obtained from the tabular values 
using a third-order Lagrange interpolation 24 stored on tape for use in 
evaluating the integral (12) on the IBM 7090. 

The phase integral S' is put into suitable form for machine computa- 
tion by dividing each of three axes, £, i? and .r/a_into equally spaced inter- 
vals. The intervals along the three axes arc A£, At?, and A(.r/a) respec- 
tively. A point in the £-77 plane now becomes a point in a grid with 
coordinate (wiA£, W2A77), where mi and m 2 are integers. The value of the 
integral (12) is now approximated by a summation over the integers 
mi and mj 

S E t iV^WTT^W) s -^# «* * @ *» (14) 

m 7T() m2 =o ^ 7 sm 7rmiA£ \ a / 

The Fourier integral (12) is thus approximated by a Fourier series 
(14). This raises the question of how well the series converges to the 
integral, which in turn is determined by the size of the intervals A£ and 
At; and by Wi and m 2 . Since a Fourier series is periodic, it not only 
approximates the integral in the range of the function »S' but produces 
repetitive images outside the range^ 

The importance of the intervals a| and Ay in the back focal plane lies 
in the fact that they set a limit to the information available at the_image 
plane. High fidelity of the image point amplitudes requires that A£ and 
A77 be as small as possible relative to the amplitude detail in the back 
focal plane. Along the 7? axis, the amplitude is a monotonic decreasing 
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function with gradual changes in slope. Along the £ axis the situation is 
quite different, as seen in Fig. 7, where IBM 7090 microfilm plots of the 
diffraction amplitude or Fourier coefficient in (lOd) are shown for 
n 1 = 5 and ??, = 17 atoms. The interval along £ is A£ = 0.01. The half 
width of the primary maxima is (A£)i = \/n x , so that a good approxi- 
mation first requires that 

A£« (Wi = - 
nx 

or 

A£ii«l. (15) 

Between the principle maxima of spacing unity there are (n x - 2) sub- 
sidiary maxima and minima. The actual spacing of these subsidiaries is 

1 - 2(Ag) t = 1_ 
rii — 2 rii 

and at least three points are required to locate a maximum, zero, and 
then the minimum. If it is arbitrarily assumed that six points are a rea- 
sonable sampling density between adjacent maxima and minima, the 
sampling interval must be 



or 



Qrti 



A* * g • 



On this basis, then, A£ = 0.01 is adequate for n x = 5. If n x = 17, 
however, then A£ ^ 1/102 or the interval A£ is just adequate. As the 
number of atoms in the chain increases furtherthe detail in the ampli- 
tude distribution soon becomes smaller than A£ = 0.01. The interval 
A£ m 0.01 is thus considered too crude for n x > 25 atoms. This is ad- 
mittedly rather arbitrary, and a more detailed analysis might discover a 
better criterion. 

Along the tj axis, on the other hand, an interval At? = 0.1 is quite ade- 
quate and reduces the number of sampling points. Under these conditions 
the number of sample points in the back focal plane for each image 
point is 

V™*l^ = 12.7 X 10 3 ; U.x = w = 3.5. (16) 

At? 7)Z 
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Fig. 7 — Diffracted amplitude distributions at the hack focal plane as ob- 
tained from the microfilm tracer on the IBM 7090 for a sampling interval A£ = 
0.01; the grating spacing is a = 2 A: (a) for a chain of 5 nickel atoms; (b) for a 
chain of 17 nickel atoms. 
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The remaining consideration relating to the sample intervals is 
aliasing. This behavior can be illustrated by returning to (14). Writing 
the coefficient as <£(£, 77), the series sums at (x/a) = to 

Sr-i'(O) = E E *(**). (17a) 

»i,=0 m 2 =0 

If the sampling interval along x is A (x/a), then the sum at an image 
point x = jA(x/a) with j an integer is 

S t J (j& (£)) = m E E *(**) cos 2tt,A f |1 m,Al. ( 17b) 

If 27r/A(.r/a)miAi = multiple of 2x then the cosine terms are unity 
and (17a) and (17b) are equal. Thus 

jA (- J A£ = integer. (17c) 

The aliasing period 2 * is then 

1 

Aj = 



© 



A - A* 



or 



alias period = =. (17d) 

Thus the critical sampling interval for these computations is A£, since 
it determines both the resolution (15) and the aliasing period. The 
effect is well illustrated in Fig. 8, comparing the profiles for n x = 5 with 
a| = 0.01 and a| = 0.1. The repeating nature of the image or aliasing 
is evident for a| = 0.1 when the aliasing period is only ten. With A£ = 
0.01, the period is 100 and not seen on the plot. If the number of atoms 
ni were now increased to nine the image profile with A£ = 0.1 would 
show no break at the end of an atom chain and so would appear to be 
an infinite chain. 

The sampling interval along the x axis in the image plane has an 
effect on the representation of the finer details in the profile. When the 
detail approaches the interval length A (x/a) the representation be- 
comes inaccurate with "bumps" and "angles" rather than a smooth 
curve. This is illustrated in Fig. 9 for the case of a single atom, ni = 1, 
for which 

= E E /(V(m,Al) 2 + (m^rj) 2 ) cos 2vjA ( - J m,Al 

,„,=U „, 2 =0 VV 



15.5 




















14.5 














(a) 






1J.5 
















































































9.5 
8.5 

7.5 
6.5 


















































*-4A-» 
























4.5 
3.5 
2.5 
1.5 


















































































IT 




















1 1 










-0.5 




















Fig. 8 — Computer-drawn image plane amplitude profiles along the axis x° 
of the nickel atom chain. For this case a = 2 A, «i_ = 5 and C„' = (u a 2 ) — in 
equation (12a). The sampling intervals in (a) are A£ = 0.01 and A(x°/a) = A77 = 
0.1. In (b) the sampling interval A£ = 0.1 and the profile shows an aliasing period 
of (AT) -1 = 10. The aliasing period in (a) is 100. 
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Fig. 9 — Profiles for equation (12) with n\ = 1 nickel atom, C 6 ' = and upper 
hounds {max = inmx = 2.5 diffraction maxima; A£ = 0.01 and Ay = 0.1. (a) Inter- 
val A(x/a) = 0.1 resulting in unresolved structure near origin, (b) Interval A(x/a) 
= 0.01 with expanded scale showing resolution of structure . 
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Fig. 10 — Effect of upper bound on AS': ni = 1 nickel atom and C 6 '_= 0. The 
upper bound is 1.5 diffraction maxima and A(x/a) = At; = 0.1 and A£ = 0.01. 
Compare with Fig. 9. 

with Af = 0.01, At; = 0.1 and A(x/a) = 0.1. The profile is shown in 
Fig. 9(a). The number of points involved in this profile is 

Wx .fm^c 20_ = gg^ x 1Q B 

At? A£ ^(x 



The neighborhood of the central maximum shows irregularities that are 
smoothed out to a more faithful representation when A(x/a) = 0.01 as 
seen in Fig. 9(b). The number of points involved in the profile is the 
same since the range has been reduced by a factor of ten. 

The subsidiary maxima and minima near the principal maximum of 
Fig. 9(b) are due to aperture diffraction and would smooth out if the 
upper limit of the integral were extended to infinity. The radius of the 
atom in the profile is not distinct, since the shoulder of the maximum will 
approach the axis asymptotically. 

The size of the aperture given by the upper bounds £ ma x and T/ mox is 
reduced from 3.5 in Fig. 9 (b) to 1.5 in Fig. 10. The reduction in aperture 
results in a loss of contrast to about 8 per cent. The upper bound 1.5 
corresponds to an objective half-angle of 2.8 X 10~ rad. 
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